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This  report  uses  "the  Doubly  Asymptotic  Approximation  to  consider  the  shock 
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In  particular,  we  estimate  the  cutoff  kinetic  energy  for  an  infinite  unrestraine 
air-backed  flat  plate  on  which  an  exponentially  decaying  acoustic  wave  is 
incident.  A  comparison  is  given  for  the  titanium  and  steel  at  equal  plate 
thickness.  Similarly,  a  comparison  is  given  for  conventional  warheads  and 
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structural  response  times.  As  expected,  there  is  also  considerable 
variation  in  the  energy  developed  in  the  plate. 
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FOREWORD 


''  ;>The  ultimate  objective  of/fbhe  work  presented- hery  i s  the  direct 
prediction  of  the  onset  and  details  of  submarine  pressure  hull  rupture  caused  by 
underwater  explosive  attack.  The  immediate  objective  of  thfrworlpis  to  bring 
out  some  of  the  implications  of  an  important  analytical  tool  known  as  the 
Doubly  Asymptotic  Approximation.-.  The  results -reported  here  provide  an  indication 
of  the  effect  of  using  light  structural  materials  such  as  titanium  and  of  using 
nonconventional  warheads  such  as  copper-lined  shaped  charge  warheads.  / 


/TJ94* 
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INTRODUCTION 


The  ultimate  objective  of  the  present  work  is  the  prediction  of  submarine 
structural  response  to  underwater  explosions.  A  major  element  in  the  physics  Is 
fluid-structure  interaction  by  virtue  of  which  pressures  exerted  by  the  fluid  are 
relieved  as  the  structure  complies.  Even  though  fluid-structure  interaction 
comprises  a  formidable  modeling  problem,  there  has  been  some  success  treating  it 
using  a  simplification  known  as  the  Doubly  Asymptotic  Approximation  (DAA)  [1]. 
Here  we  illustrate  some  consequences  of  the  DAA  in  a  model  problem  involving  an 
exponentially  decaying  acoustic  wave  impinging  on  an  infinite  air-backed  unre¬ 
strained  flat  plate,  as  illustrated  in  Figure  1. 

DOUBLY  ASYMPTOTIC  APPROXIMATION 

For  the  model  problem  the  incident  pressure  is  subject  to  exponential  decay 
typical  of  underwater  explosions: 

Pj  =  pQ  exp(-t/e) 

The  total  pressure  on  the  plate  may  be  thought  of  as  the  sum  of  the  incident 
pressure,  the  reflected  pressure,  and  the  radiated  pressure  arising  from  plate 
compliance: 

PT  =  Pi  +  Pr  +  Pr 
The  DAA  uses  the  approximations 

"r  •  Pi  Pr  *  -»c* 

where  p  is  the  water  density,  c  is  the  acoustic  wave  speed,  and  v  is  the  plate 
velocity. 

The  equilibrium  equation  becomes 
dv 

m3t  *  PT 

*  2p0  exp(-t/e)-Pcv 

*  Geers,  T.  L.,  "Transient  Response  Analysis  or  Submerged  Structures,"  in 
Finite  Element  Analysis  of  Transient  Nonlinear  Structural  Behavior, 

A$m£,  1975. - -  - 

2 

Kennard,  E.  H.,  "The  Effect  of  Pressure  Wave  on  a  Plate  or  Diaphragm,"  in 
Underwater  Explosion  Research.  Office  of  Naval  Research,  1950. 

3 

Taylor,  G.  I.,  "The  Distortion  Under  Pressure  of  a  Diaphragm  which  is  Clamped 
Along  its  Edge  and  Stressed  Beyond  its  Elastic  Limit,"  in  Underwater  Explosion 
Research,  Office  of  Naval  Research,  1950. 
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where  m  is  the  plate  mass  per  unit  area.  Now  using  the  dimensionless  quantities 

w  ■  v/c  t  =  t/e 

k  *  pce/m  v  *  Zpoe/cm 

the  equation  (1)  may  be  rewritten 

^  +  m  »  v  exp(-t).  (2) 

With  the  initial  condition  w(o)  a  o,  the  solution  of  (2)  is 

w  •£-  [exp(-KT)-exp(-c)]  (3) 

1  * 

KINETIC  ENERGY  AT  CUTOFF 

Cutoff  occurs  when  the  total  pressure  due  to  the  shock  loading  vanishes.  We 
assume  that  no  shockwave  energy  is  transferred  to  the  plate  after  cutoff.  However, 
in  reality,  a  considerable  amount  of  energy  may  be  transferred  by  afterflow,  bubble 
pulses,  etc.  This  has  to  be  the  subject  of  a  subsequent  analysis.  The  kinetic 
energy  at  cutoff  should  be  considered  the  shock  wave  energy  available  for  damage. 
Pressure  cutoff  therefore  occurs  when 

pcv  *  2p0  exp(-t/e) 

and  hence  when 

kw  *  v  exp(-t).  (4) 

Substituting  (4)  into  (3)  and  solving  for  the  cutoff  time  tc  leads  to 


rc  =  -{en  tcj'(i-ic) 
The  cutoff  velocity  is  w£  where 


w c  *  7  exp(-t  ) 

Most  importantly,  the  kinetic  energy  at  cutoff  is 


K  * i# V 2 
C  2  C 


*  v  e  ic  * mK  (5) 

where  u  *  2p2Q/pc 

The  prominence  of  the  quantity  <  *  cce/m  illustrates  the  physics  very  succinctly 
It  is  nothing  but  the  ratio  of  pressure  decay  times  to  structural  (inertial) 
response  times. 
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COMPUTATIONAL  RESULTS 


Computations  are  now  presented  for  two  cases. 

Case  1:  Here  the  time  constant  0  Is  held  fixed  and  only  the  mass  m  is  varied.  This 
corresponds  to  comparisons  of  the  effect  of  using  titanium  in  place  of  steel. 

An  appropriate  dimensionless  counterpart  of  the  cutoff  kinetic  energy  is  given 
by  xi  *  Kc/u0,  for  which 

i+ie 

X,  -  <  ** 

The  limiting  cases  of  are 


lim  xx  =  0  lim  xt  =  0 

tc-K) 

lim  Xj  *  e”2  =  .135 

K-i 

The  maximum  value  of  xt  corresponds  to  k  =  1.  It  should  be  recalled  that 
k  -  pce/m  and  so  <  is  inversely  proportional  to  m. 

The  dependence  of  kinetic  energy  at  cutoff  on  plate  mass,  represented  by  the 
function  x^k),  is  shown  in  Figure  2.  Illustrative  examples  of  these  results  will 
be  given  in  the  next  section. 

Case  2:  Here  the  mass  m  is  fixed,  while  only  the  time  constant  e  is  allowed  to 
vary.  Now  a  dimensionless  counterpart  of  the  cutoff  kinetic  energy  is  x2  =  K  /X 
where  x  *  2p2Qm/p2c2,  for  which  c 

x2  =  K  2/1 -k 

This  quantity  has  the  limiting  cases 

11m  x2  *  0  lim  x2  *  1 


lim  x2  *  e'2  =  .135 
k-m 

Unlike  xx»  the  maximum  of  x2  occurs  at  k-k»  and  not  at  *  *1. 

The  dependence  of  kinetic  energy  at  cutoff  on  decay  time  e,  represented  by  the 
function  X2(*)*  Is  shown  in  Figure  3.  Recall  that  <  is  directly  proportional  to  9. 
Illustrative  examples  of  these  results  will  be  given  in  the  next  section. 
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EXAMPLES 


For  the  conventional  UNDEX  explosive  represented  by  50  pounds  of  TNT  at  300 
feet  standoff,  reference  2  gives  the  value 

•4 

*  7*10  sec- 

For  a  copper-lined  shaped  charge  warhead  fired  under  water,  we  found  some  data 
from  which  we  estimated  the  decay  time  as 

-5 

0,  *  3.5-10  sec- 
2 

For  a  one  inch  thick  steel  plate 

(m/pc)^  *  1.3-10’^  sec 


while  for  a  one  inch  thick  titanium  plate 

(m/pc)2  =  0.78*10"Sec 
We  may  form  the  following  table: 


Conventional 
K »Xi  /• 135 ,  X2 

steel  5.38,  0.638,  0.464 

titanium  8.97,  0.476,  0.579 


Copper  Jet 
k  «x 1  /  • 135 ,  X2 

0.27,  0.759,  0.028 
0.45,  0.903,  0.055 


The  values  of  x^-135  suggest  that  a  titanium  flat  plate  develops  about  15-30% 
less  cutoff  kinetic  energy  than  a  steel  flat  plate  with  the  same  plating  thicknesses. 
The  values  of  x2  suggest  that,  for  both  metals  and  a  fixed  incident  pressure 
amplitude  p  ,  tne  cutoff  kinetic  energy  for  a  shaped  charge  jet  is  a  small  fraction 
of  that  for  a  conventional  warhead. 


DISCUSSION 


The  Doubly  Asymptotic  Approximation  has  been  used  to  investigate  the  consequences 
of  certain  target  and  warhead  variations  on  the  shock  wave  energy  available  for 
damage.  At  similar  plate  thickness,  a  titanium  structure  may  develop  15-30%  less 
energy  than  a  steel  structure.  At  the  same  incident  pressure  amplitudes,  a 
copper-lined  shaped  charge  warhead  may  introduce  only  a  small  fraction  of  the  energy 
introduced  by  a  conventional  warhead. 
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FIGURE  3  EFFECT  OF  TIME  CONSTANT  ON  THE  KINETIC  ENERGY  AT  CUTOFF 
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